Winds from accretion disks driven by the radiation and magnetocentrifugal force. 
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ABSTRACT 



We study the two-dimensional, time-dependent hydrodynamics of radiation-driven winds from luminous 
accretion disks threaded by a strong, large-scale, ordered magnetic field. The radiation force is mediated 
primarily by spectral lines and is calculated using a generalized multidimensional formulation of the Sobolev 
approximation. The effects of the magnetic field are approximated by adding into the equation of motion 
a force that emulates a magnetocentrifugal force. Our approach allows us to calculate disk winds when the 
magnetic field controls the geometry of the flow, forces the flow to corotate with the disk, or both. In particular, 
we calculate models where the lines of the poloidal component of the field are straight and inclined to the disk 
at a fixed angle. 

Our numerical calculations show that flows which conserve specific angular velocity have a larger mass 
loss rate than their counterparts which conserve specific angular momentum. The difference in the mass 
loss rate between these two types of winds can be several orders of magnitude for low disk luminosities but 
vanishes for high disk luminosities. Winds which conserve angular velocity have much higher velocities than 
angular momentum conserving winds. Fixing the wind geometry stabilizes winds which are unsteady when 
the geometry is derived sclf-consistently. The inclination angle between the poloidal velocity and the normal 
to the disk midplane is important. Non-zero inclination angles allow the magnetocentrifugal force to increase 
the mass loss rate for low luminosities, and increase the wind velocity for all luminosities. The presence of the 
azimuthal force does not change the mass loss rate when the geometry of the fiow is fixed. 

Our calculations also show that the radiation force can launch winds from magnetized disks. The line force 
can be essential in producing magnctohydrodynamical (MHD) winds from disks where the thermal energy is 
too low to launch winds or where the field lines make an angle of < 30° with respect to the normal to the 
disk midplane. In the latter case the wind will be less decoUimated near its base by the centrifugal force 
and its coUimation far away from the disk can be larger than the coUimation of its centrifugally-driven MHD 
counterpart. 



1. Introduction 

Accretion disks are believed to lose mass via powerful outflows in many astrophysical environments such as active 

galactic nuclei (AGN); many types of interacting binary stars, non-magnetic cataclysmic variables (nMCVs), for instance; 
and young stellar objects (YSOs). Magnetically driven winds from disks are the favored explanation for the outflows in 
many of these environments. Blandford & Payne 1982 (see also Pelletier & Pudritz 1992) showed that the centrifugal 
force can drive a wind from the disk if the poloidal component of the magnetic field. Bp makes an angle of > 30° with 
respect to the normal to the disk surface. Generally, centrifugally-drivcn MHD disk winds (magnetocentrifugal winds for 
short) require the presence of a sufficiently strong, large-scale, ordered magnetic field threading the disk and the poloidal 
magnetic field to be comparable to the toroidal magnetic field, \B0/Bp\ ^ 1 (e.g., Cannizzo & Pudritz 1988, Pelletier & 
Pudritz 1992). Additionally, magnetocentrifugal winds require some thermal assistance to flow freely and steadily from 
the surface of the disk, to pass through a slow magnetosonic surface (e.g., Blandford & Payne 1982). 

Many authors have studied magnetocentrifugal winds from a Keplerian disk (e.g., Ouyed & Pudritz 1997, Ustyugova 
et al. 1999, Krasnopolsky, Li & Blandford 1999, and references therein). These studies are either analytic, looking for 
stationary, often self-similar solutions or numerical, looking for both stationary and time-dependent solutions. However 
in numerical simulations of magnetocentrifugal winds, the lower boundary is between the slow magnetosonic surface and 
the Alfven surface. This specification requires setting the mass loss rate in advance (e.g., Bogovalov 1997). Simulations of 
magnetocentrifugal winds considering the whole disk (even regions below the slow magnetosonic surface) and not requiring 
an ad hoc mass loss rate, are only now becoming feasible. This requires, however, an accurate treatment of the radiation 
and gas pressure effects, among other physical processes. 

Thermal expansion and the radiation force have been suggested as other mechanisms that can drive disk winds. These 
mechanisms can produce powerful winds without the presence of a magnetic field. Winds arc likely thermally driven in 
X-ray-irradiated accretion disks in systems such as X-ray binaries and AGNs (e.g., Bcgelman, McKee & Shields 1983, 
Woods et al. 1996). These winds require the gas temperature to be so high that the gas is not gravitationally bound. In 
such cases, the radiation driving is probably not important because the gas is fully ionized, at least in the hottest regions, 
and the radiation force is only due to electron scattering. However in the regions where gas is cooler and not fully ionized 



2 



the radiation force will be enhanced by spectral lines and play an important role in controlling the dynamics of the flow. 
In fact Murray et al. (1995) designed a line-driven disk wind model specifically for AGNs. 

Radiation-driven disk winds have been extensively modeled (e.g., Pereyra, Kallman & Blondin 1997; Proga, Stone & 
Drew 1998, hereafter PSD I; Proga 1999; Proga, Stone & Drew 1999, hereafter PSD II). These recent studies showed that 
radiation pressure due to spectral lines can drive winds from luminous disks. This result has been expected (e.g., Vitello 
& Shlosman 1988). These studies, in particular those by PSD I, also showed some unexpected results. For example, the 
flow is unsteady in cases where the disk luminosity dominates the driving radiation field. Despite the complex structure of 
the unsteady disk wind, the time-averaged mass loss rate and terminal velocity scale with luminosity, as do steady flows 
obtained where the radiation is dominated by the central object. 

To calculate the line force for disk winds, PSD I adopted the method introduced by Castor, Abbott & Klein (1975, 
hereafter CAK) and further developed by Friend & Abbott (1986), Pauldrach, Puis & Kudritzki (1986) for one-dimensional 
radial flows within the context of a wind from hot, luminous OB stars. To extend the CAK method to model multi- 
dimensional disk winds, it is necessary to accommodate the effects of the three-dimensional velocity field and the direction- 
dependent intensity. Owocki, Cranmer & Gayley (1996) showed that these effects can lead to qualitatively different results 
compared to those obtained from a one-dimensional treatment in the case of a rapidly rotating star. The most difficult 
aspect of calculating the line force due to a disk is in the evaluation of the integral involving the velocity gradient 
tensor over the entire solid angle occupied by the radiating surface. PSD I used an angle-adaptive quadrature to ensure 
an accurate result. However, computational limitations required that they simplified the integrand, retaining only the 
dominant terms in the velocity gradient tensor. 

PSD II generalized the PSD I method and introduced a new quadrature that avoids any simplification of the integrand. 
This allowed them to evaluate the radiation force for completely arbitrary velocity fields within the context of the CAK 
formalism. They applied this method to recalculate several disk wind models first discussed in PSD I. These more 
physically accurate models show that PSD I's more approximate method was very robust. The PSD II calculations 
predict total mass loss rates and velocities marginally different from those published in PSD I. Additionally, PSD II find 
that models which display unsteady behavior in PSD I are also unsteady with the new method. The largest change caused 
by the new method is in the disk-wind opening angle: winds driven only by the disk radiation are more polar with the 
new method while winds driven by the disk and central object radiation are typically more equatorial. 

In this paper, we will numerically check how strong, ordered magnetic fields can change disk winds driven by the line 
force for a given disk luminosity. We assume that the transport of angular momentum in the disk is dominated by local 
disk viscosity, for instance due to the local shear instability in weakly magnetized disk (Balbus & Hawley 1997). Instead 
of adding the magnetic fields to the PSD II model and solving consistently the equations of MHD, we simply start by 
adding some of the effects due to the magnetic fields, namely (1) the azimuthal force so the wind conserves the specific 
angular velocity along the streamlines and (2) the force perpendicular to the field lines so the geometry of the disk wind 
is controlled by the magnetic field. We thus adopt the popular concept that the magnetic field dominates outside the 
disk and at least near the disk surface, one can think of the magnetic field lines as rigid wires that control any flow (e.g., 
Blandford & Payne 1982; Pelletier & Pudritz 1992). Such an approach is clearly simplistic and can be applied only to 
sub-Alfvcnic flows. Nevertheless, the results presented here provide a useful exploratory study of line-driven disk winds 
in the presence of a strong, large-scale field magnetic threading the disk. 

The organization of the paper is as follows. In Section 2 we describe our numerical calculations; in Section 3 we present 
our results; in Section 4 we conclude with a brief discussion. 



2. Method 

To calculate the structure and evolution of a wind from a disk, we solve the equations of hydrodynamics 

^+PV.V = 0, (1) 

P^=-V(pa+pg + pF™'^ + pF'"'= (2) 

where p is the mass density, v the velocity, g the gravitational acceleration of the central object, F''"'^ the total radiation 
force per unit mass, and F™"^ the total 'magnetocentrifugal' force. The term with F™"^ ensures that gas flows along an 
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assumed direction and conserves its specific angular velocity. The gas in the wind is isothermal with a sound speed Cg. 
We solve these equations in spherical polar coordinates (p, 9, (j)). 

The geometry and assumptions needed to compute the radiation field from the disk and central object are as in PSD II 
(see also PSD I). The disk is flat, Kcplcrian, geometrically-thin and optically-thick. Wc specify the radiation field of the 
disk by assuming that the temperature follows the radial profile of the so-called a-disk (Shakura & Sunyaev 1973), and 
therefore depends only on the mass accretion rate in the disk, M^, and the mass and radius of the central object, M* 
and r*. In particular, the disk luminosity, = GM^Ma/2r,:. In models where the central object radiates, wc take into 
account the stellar irradiation of the disk, assuming that the disk re-emits all absorbed energy locally and isotropically. 
We express the central object luminosity in units of the disk luminosity = xLjj- See PSD I and PSD II for further 
details. 

Our numerical algorithm for evaluating the line force is described in PSD II. For a rotating flow, there may be an 
azimuthal component to the line force even in axisymmetry. However we set this component of the line force to zero 
because it is rather weak as compared to other components and is not of great importance (e.g., PSD II). Then wc assume 
that the rotational velocity is determined by the azimuthal component of the magnetocentrifugal force. The description 
of our calculation of the 'magnetocentrifugal' force follows. 

Wc choose the simplest geometry of the magnetic field and the flow: the poloidal component of the magnetic field. 
Bp and of the velocity, Vp are parallel to one another, and the inclination angle i between Vp and the disk midplane is 
fixed for all locations and times. In other words, we constrain the geometry of the flow to straight cones of pre-specified, 
radius- and time-independent opening angle. In practice, wc impose this geometry in the following way: (i) we evaluate 
the physical accelerations and the curvature terms in eq. 2, hereafter collectively referred to as 

F' = -^V(pc2)+g + F™'^ + F\ (3) 

where represents the curvature terms appearing on the left hand side of eq. 2 when is expressed in the spherical 
polar coordinate system (e.g., Shu 1992): 



V 

(ii) we calculate the component of F' perpendicular to a streamline: 



(4) 



F'± = AF', (5) 

where 

/ sin^ {9 - i) sin (6* - i) cos [9 - i) \ 

A = sin (0 - i) cos {9 - i) cos^ {9 - i) , (6) 

\ / 

(iii) and finally we subtract F^ from F' so the total remaining force is tangent to the streamline. 

The disk outflow that conserves its specific angular velocity, Q. = v^/{rsm6) along a streamline is simply the flow that 
corotates with the disk, i.e., is constant and equals to the disk rotational velocity, = w</,(rD, 90°)/rD at the footpoint 
of the streamline on the disk at the radius, . To ensure corotation of the wind, we introduce the azimuthal force: 

F* = I 0_ 1 . (7) 



Summarizing, our magnetocentrifugal force is 



2 Vej, (Vr sin 9-\-V0 COS 0) 

r sin 



F'"'^ = -AF'-I-F*. (8) 



As in PSD II, we use the ZEUS-2D code (Stone & Norman 1992) to numerically integrate the hydro dynamical equations 
(1) and (2). We did not change the transport step in the code and so the changes we made are just to the source step as 
outlined above. 
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Wc explore disk wind models for three cases: (I) a wind corotates with the disk (F'^ = and F^* = 2v^{vr sinO + 
ve cos 9) /{r sin 9)) ; (II) the wind geometry is fixed by adopting a fixed inclination angle between the poloidal component 
of the velocity and the disk midplane ( F'_l = AF' and F* = 0) ; and (III) a combination of (I) and (II) ( F'_l = AF' and 
F0* = 2v^{vr sin 9 + vg cos 9) / {r sin 9)). 

In case I with a; = 0, we decided to put an additional constraint on the radial velocity to produce outward flows. Near 
the disk midplane, fluctuations of the density and velocity can occur (e.g., PSD I). Additionally, the radial component 
of the radiation force is negative near the disk midplane and the central object (e.g., Icke 1980, PSD I). Thus Vr can 
be less than zero; gravity can dominate and pull the flow toward the center. To avoid this we apply the condition: 
Vr = max{vr, — cot9ve). This condition does not change the mass loss rate, but does reduce the wind opening angle. 

To produce transonic flows in cases II and III, we need to increase the density along the disk midplane, po, from 
10~^ g cm~"^ (as used by PSD I) up to 10~"^ g cm~^. For low po, the wind velocity does not depend on p, the mass loss 
rate is proportional to po and may be higher than the mass accretion rate. For po ^ 10~* g cm~^ the mass loss rate is 
dramatically lower, does not change with po, and the flow is subsonic near the disk midplane. These are desired properties 
of the outflow; we expect gas near the midplane to be nearly in hydrostatic equilibrium. Equally important, we assume 
that the disk is in a steady state so the mass loss rate should be lower than the mass accretion rate. 

In reality, gas near the disk midplane (inside the disk) is in hydrostatic equilibrium for the densities lower than those 
we assumed because the radiation force is lower as the radiation becomes isotropic. We would like to stress that we treat 
the region very close to the midplane as a boundary condition and do not attempt to model the disk interior. 

As in PSD I and PSD II, we calculate disk winds with model parameters suitable for a typical nMCV (see PSD I's 
Table 1 and our Table 1). We vary the disk and central object luminosity and the inclination angle. We hold all other 
parameters flxed, in particular the parameters of the CAK force multiplier: k = 0.2, a = 0.6 and M^ax = 4400 (see 
PSD I). Nevertheless we can use our results to predict the wind properties for other parameters and systems such as 
AGN - by applying the dimensionless form of the hydrodynamic equations and the scaling formulae as discussed in PSD I 
and Proga (1999). 

3. Results 

PSD I and PSD II showed that radiation-driven winds from a disk fall into two categories: 1) intrinsically unsteady 
with large fluctuations in density and velocity, and 2) steady with smooth density and velocity. The type of the flow is set 
by the geometry of the radiation fleld, parametrized by x: if the radiation fleld is dominated by the disk {x < 1) then the 
flow is unsteady, and if the radiation is dominated by the central object {x ^ 1) then the flow is steady. The geometry 
of the radiation field also determines the geometry of the flow; the wind becomes more polar as x decreases. However 
the mass-loss rate and terminal velocity are insensitive to geometry and depend more on the total system luminosity, 
Ld + L^. We recalculated some of the PSD II models to check how inclusion of the magnetocentrifugal force will change 
line- driven disk winds. 

Figure 1 compares the density in the wind in two models from PSD II where they used a generalized CAK method (top 
panels), models using PSD IPs method but conserving specific angular velocity - our case I (middle panels), and models 
using PSD IPs method, conserving specific angular velocity and having fixed inclination angle of the stream lines - our 
case III (bottom panels). Models in case III (bottom panels) are with the inclination angles approximately equal to those 
in the PSD II case (the top panels). The left column shows the results for a model with Ma = 10~* yr~^ and a; = 0, 
while the right column shows the results with Ma = tt x 10~^ yr"-'^ and x = \. The former corresponds to the fiducial 
a; = model while the latter corresponds to the fiducial a; = 1 model discussed in detail in PSD I and PSD II. 

We start with describing the a; = wind model (left column panels). The most striking yet expected difference between 
PSD IPs case (top) and our case I (middle) is a decrease of the wind opening angle, w from 50° to 15° when the specific 
angular momentum conservation is replaced by the specific angular velocity conservation. This result is caused by a much 
stronger centrifugal force in our case I than in the PSD II case which is also reflected in an increase of the radial velocity 
by more than one order of magnitude (Table 1). Another big difference between the two cases is in the mass loss rate, 
which increases from 5.5 x 10~^* Mq yr~^ to 1.3 x 10~^^ yr~^. An unchanged property in these two cases is that 
the wind is unsteady. The quantitative changes between PSD I's case (top) and case III (bottom) are similar or smaller 
than the changes described above (see Table 1). The most pronounced change caused by holding the inclination angle 
fixed is that the wind becomes steady. However this is not surprising because models with the fixed, rigid geometry are 
pseudo one-dimensional and the flow has one degree of freedom in space. In PSD I and PSD II models where the wind 
geometry is calculated self-consistently, a strong radial radiation force is required to produce a steady outflow. 
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The x=l wind model (top right hand panel of Figm'e 1) is an example of a steady outflow from PSD II. This model 
remains steady in the two other cases shown in the figure and in case I listed in Table 1 (run Ic). The model changes 
from case to case mainly in the radial and rotational velocities which increase when the azimuthal force is added (middle 
and bottom panel). For example, the radial velocity increased from 3500 km s^-'^ in PSD IPs run C (top panel) to 28000 
km s~^ in our run IIIc (bottom panel). As in the x=0 wind model, adding the azimuthal force also reduces the wind 
opening angle (middle panel). However an increase in the mass loss rate from PSD IPs case (top panel) to our case II 
(bottom panel) is only a factor of 2. The mass loss rate is practically the same in case II and case III (run IIc and IIIc). 

Models presented in Figure 1 both for a; = and x = 1 suggest that an increase in the mass loss rate due to conservation 
of the specific angular velocity decreases with the disk luminosity. Our other results, shown in Table 1, confirm this. For 
example, the x=0 model lb with Ma = tt x 10~^ Mq yr"-'^ has higher than the corresponding PSD IPs model B by a 
factor of ~ 3 while an increase of Myj for models with = 10~^ Mq yr~^ is by a factor ~ 24 (i.e., runs A and la). 

The mass loss rate of the line-driven disk wind also increases by imposing the wind geometry such that i > 0" as 
comparison of our results in case II and from PSD II reveals. Once the geometry is fixed the mass loss rate does not 
change with adding the azimuthal force - the corresponding models in our case II and III have the same Mw- 

Purely line-driven x=0 winds have streamlines perpendicular to the disk midplane near the midplane because the total 
horizontal force is negligible in comparison with the vertical force due to lines. This property of disk winds has been 
assumed in analytic studies of line-driven disk winds (e.g., Vitello & Shlosman 1988). We calculated a model assuming 
i = 0° (run Ila) and find that the mass loss rate is very similar to the corresponding PSD IPs run A. This confirms 
then that the mass loss rate in line-driven disk winds for a; = is determined in the regions where the matter flows still 
perpendicularly to the disk. 



4. Conclusions and Discussion 

We have studied winds from accretion disks with very strong organized magnetic fields and the radiation force due to 
lines. We use numerical methods to solve the two-dimensional, time-dependent equations of hydrodynamics. We have ac- 
counted for the radiation force using a generalized multidimensional formulation of the Sobolev approximation. To include 
the effects of the strong magnetic field we have added to the hydrodynamic equation an approximate magnetocentrifugal 
force. Our approach of treating the latter allows us to study cases where the wind geometry is controlled entirely by the 
magnetic field with the lines of the poloidal component of the field. Bp being straight and inclined at a fixed angle to the 
disk midplane where the specific angular velocity is conserved along the streamlines due to the magnetic azimuthal force, 
or both. 

Our simulations of line-driven disk winds show that inclusion of the azimuthal force, F* which makes the wind corotate 
with the disk, increases the wind mass loss rate significantly only for low disk luminosities. As expected such winds have 
much higher velocities than winds with zero azimuthal force. Fixing the wind geometry effectively reduces the flow from 
two-dimensional to one-dimensional. This in turn stabilizes winds which are unsteady when the geometry is allowed to 
be determined self-consistently. The inclination angle between the poloidal velocity and the normal to the disk midplane 
is important. If it is higher than zero it can significantly increase the mass loss rate for low luminosities, and increase the 
wind velocity for all luminosities. The presence of the azimuthal force does not change the mass loss rate if the geometry 
is fixed. 

It is very intriguing that the mass loss rate can be enhanced by the magnetocentrifugal force the most for low luminosities 
where line force is weak and at the same time the mass loss rate remains a strong function of the disk luminosity for all 
luminosities. In particular, there is no wind when Ma ^ 10~* M0yr~^ as in purely line-driven case (e.g., PSD I; Proga 
1999). 

Model Ila' with i = is steady and has the same Mu, as the corresponding PSD II model. The mass loss rate is 
then the same regardless of the time behavior. This is consistent with PSD I conclusion that the mass loss rate depends 
predominately on the total system luminosity. 

Drew & Proga (1999) applied results from PSD I, PSD II and Proga 1999 to nMCV. In particular, they compared 
mass loss rates predicted by the models with observational constraints. They concluded that either mass accretion rates 
in high-state nMCV are higher than presently thought by a factor of 2-3 or that radiation pressure alone is not quite 
sufficient to drive the observed hypersonic flows. If the latter were true then an obvious candidate to assist radiation 
pressure in these cases is MHD (e.g.. Drew & Proga 1999). 

An increases of the mass loss rate due to inclusion of a magnetocentrifugal force (our case I) brings our predictions close 
to the observational estimates for nMCVs. However we should bear in mind that at the same time, the radial velocity of 
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The increase 

of the wind velocity may not be large if the magnetic field is moderately strong and spins up the disk wind only close to 
its base where the mass loss rate is determined. Then farther away from the disk the wind would gain less or no angular 
momentum and velocities would be similar to those in the case without the magnetic field at all. 

Outflows which conserve specific angular velocity have a higher rotational velocity than those which conserve specific 
angular momentum. The former also have higher terminal velocities due to stronger centrifugal force. Additionally in 
the wind where the specific angular velocity is conserved the rotational velocity is comparable to the terminal velocity 
while in the wind where the specific angular momentum is conserved the rotational velocity decreases asymptotically to 
zero with increasing radius and is therefore much lower than the terminal velocity. Thus we should be able to distinguish 
these two kinds of winds based on their line profiles. For example, highly rotating winds should produce emission lines 
much broader than slowly rotating, expanding winds, if we see the disk edge-on. 

So far we have discussed changes of the line-driven wind caused by inclusion of the magneto centrifugal force. However 
we can also anticipate some changes in the magnetocentrifugal wind if we add the line force. For example, there may be 
a difference in the wind collimation caused by the line force. The centrifugal force decollimates the magnetocentrifugal 
outflow near the disk as i must be greater than 30°. In such a case the magnetic field coUimates the outlow only beyond 
the Alfvcn surface where the pinch force exerted by the toroidal component of the field operates. In the case with a 
significant line driving, the magnetic field does not have to be inclined at the angle > 30°. Then in such case, outfiows 
are more collimated than pure magnetocentrifugal outflows from the start and may end up more collimated far away from 
the disk. 

Let us now discuss some results of other related works. Recently Ogilvie & Livio (1998) studied a thin accretion disk 
threaded by a poloidal magnetic field. Their purpose was to determine qualitatively how much thermal assistance is 
required for the fiow to pass through the slow magnetosonic surface. They found that a certain potential difference must 
be overcome even when i > 30° and that thermal energy is not sufficient to launch an outlow from a magnetized disk. 
Ogilvie & Livio suggested that an additional source of energy, such as coronal heating may be required. PSD I showed the 
radiative line force can drive disk winds. Thus the 'missing' energy may be in the radiation field. Our calculations here 
show that if the disk luminosity is sufficiently high the line and magnetocentrifugal forces produce strong transonic winds 
regardless of the wind geometry. Then the magnetocentrifugal force can assist the line force in producing disk winds - 
our approach to the problem or the line force or/and the thermal force can well assist MHD - Ogilvie & Livio's approach 
(see also Wardle & Kongil 1993; Cao & Spruit 1994, for instance). In seeking a steady state solution of a line-driven 
disk wind, Vitello & Shlosman (1988) also found that the flow must overcome a potential difference - an increase of the 
vertical gravity component with height. They suggested that for the radiation force to increase with height above the disk 
midplane and overcome this potential difference a very particular variation in the ionization state of the gas is required. 

There are a number of limitations of our treatment of magnetic fields which are worthy of mention. 

Our models do not include the whole richness of MHD because we approximate the Lorentz force by the 'magnetocen- 
trifugal' force (eqs 3-8) instead of solving self-consistently the equation of motion and the induction equation. We have 
included effects of a magnetic field on the disk winds in a simplistic manner that mimics some effects of a very strong, 
organized magnetic field where ^ Bp near the disk surface. The; magnetic field is treated as a rigid wire that controls 
the flow geometry outside the disk. More quantitatively, this corresponds to the situation where outside the disk, at 
least near the disk photosphere, the magnetic field dominates i.e., the magnetic pressure exceeds the disk gas pressure 
f3 = 8Trp]j/Bf) < 1. Inside the disk however the situation may be different. We assume that the disk is in a steady state 
that is also stable. In such a case, the regions near the disk midplane are likely supported by the gas pressure and the 
magnetic pressure should be smaller than the gas pressure, i.e., /3 ^ 1, otherwise the disk may be unstable (Stella & 
Rosner 1984). Using the system parameters adopted here: Cs = 14 km s~^ and po = 10~^ g cm~^, we find that the /? <; 1 
condition yields a maximum value for the disk magnetic fields strentgh of the order of ~ 10^ G. Our assumption that 
poloidal magnetic field lines are straight excludes any magnetic tension in the (r, 6) plane. 

In our cases I and III we assume that the specific angular velocity, Q is conserved along a streamline from the footpoint of 
the line to infinity or rather to the outer boundary of the computational domain. However in a steady state axisymmetric 
MHD flow, the quantity which is conserved is the total angular momentum per unit mass which can be written as 

I = Ur^ sm^ 9 = flori , (9) 

where ta is the position of the the Alfven point of the flow on the streamline, and k = pVp/Bp is the mass load- another 
constant along the streamline (e.g., Pelletier & Pudritz 1992). The total angular momentum has contributions from both 
the flowing rotating gas and the twisted magnetic fleld. This means that as the material angular momentum, Or^ sin^ 6 
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increases chie to corotation with the disk, the toroidal component of the magnetic field must inc;rease so the total angular 
momentum is conserved. Then our approximation in cases I and III corresponds to the situation where near the disk 
the total angular momentum is dominated by the contribution from the twisted magnetic field and the Alfvenic surface 
is beyond our computational domain. In other words our approximation is valid in the region where the outflow is sub- 
Alfvenic. Consequently, our approach does not allow us to study coUimation of the wind by the magnetic field because 
this happens beyond the fast magnetosonic surface where the flow is super- Alfvenic. 

PSD II's models and ours in case II correspond to the other extreme case where the total angular momentum has no 
contribution from the twisted magnetic field and the material angular momentum is conserved. 

Our treatment of magnetic fields does not include effects of the magnetic pressure. The magnetocentrifugal driving 
presupposes the existence of a strong poloidal magnetic field comparable to the toroidal magnetic field near the disk 
surface, \B^/Bp\ ^ 1. However when the poloidal magnetic field is weaker than the toroidal magnetic field, \B(i,/Bp\ » 1 
the magnetic pressure may be dynamically important in driving disk winds (e.g., Uchida & Shibata 1985; Pudritz & 
Norman 1986, Shibata & Uchida 1986, Contopoulos 1995; Kudoh & Shibata 1997; Ouycd & Pudritz 1997 and references 
therein). For \B^/Bp\ » 1, there is initially the buildup of the toroidal magnetic field by the differential rotation of 
the disk that in turn generates the magnetic pressure of the toroidal field. The magnetic pressure then gives rise to a 
self-starting wind. To produce a steady outflow driven by the magnetic pressure a steady supply of the advected toroidal 
magnetic flux at the wind base is needed, otherwise the outflow is likely a transient (e.g., Konigl 1993, Contopoulos 1995, 
Ouyed & Pudritz 1997). However it is not clear whether the differential rotation of the disk can produce such a supply 
of the toroidal magnetic flux to match the escape of magnetic flux in the wind and even if it does whether such a system 
will be stable (e.g., Contopoulos 1995, Ouyed & Pudritz 1997 and references therein). 

Concluding, we would like to stress that a further development of models of radiation-driven winds from disks should 
include taking into account magnetic field but it is equally important to consider adding the radiation force to models of 
MHD winds from luminous disks. Both kinds of models are quite well understood now and if merged they could allow us 
to study better disk winds in systems such as cataclysmic variables and AGNs. As we mentioned above thermal assistance 
may not be sufficient to launch wind from a magnetized disk and we should consider not only coronal heating but also 
line driving. 
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Table 1. Summary of results for disc winds 



run 


Ma 

(Mo yr-1 ) 


X 


My, 

(Mo yr-i) 


Wr(10r*) 

(km s-i) 


CO or (90° - i)* 
degrees 


PSD II 
A 
B 
C 


10-8 

TT X 10-8 
TT X 10-8 





1 


5.5 X 
A n V 

2.1 X 


10-14 

1fl-12 

10-11 


900 
3500 
3500 


50 
60 
32 


I 

a 
b 
c 


10-8 
TT X 10-8 
TT X 10-8 





1 


1.3 X 

1.3 X 

2.4 X 


10-12 

10-11 
10-11 


15000 
20000 
32000 


15 
38 
12 


II 

a 
a' 
c 


10-8 
10-8 
TT X 10-8 






1 


6.3 X 

6.3 X 
4.2 X 


10-14 
10-13 

10-11 


600 

1100 
5000 


90t 

60t 
30t 


III 

a 
c 


10-8 
TT X 10-8 




1 


6.3 X 
4.2 X 


10-13 

10-11 


16000 
28000 


60t 
30t 



* For all models from PSD II and our models in case I, the last column lists the wind opening angle, u). While for models 

in cases II and III, the last cohimn hsts 90° — i (marked with f), where i is the assumed inchnation angle between the 
poloidal component of the velocity and the normal to the disk midplane; note that then w w (90° — i). 
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2 4 6 8 10 
r (r.) 




024-6810 024 6 810 

r (r.) r (r.) 

Fig. 1. — The color density maps in two models computed using the method of PSD 11 (top panels), the PSD II method 
with the specific angular velocity conservation (middle panels) , and the latter with a fixed inclination angle of streamlines 
(bottom panels). The left column shows the results for the x = model while the right column shows results for the 
X — 1 model (see table 2, and section 3.1 for discussion). 
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